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We show that the inclusion of an axion-like effective potential in the construction of a self- 
gravitating system of scalar fields decreases its compactness when the value of the self-interaction 
coupling constant is increased. By including the current values for the axion mass m and decay 
constant f a , we have computed the mass and the radius for self-gravitating systems made of axion 
particles. It is found that such objects will have asteroid-size masses and radii of a few meters, thus 
a self-gravitating system made of axions could play the role of scalar mini-MACHOs and mimic a 
cold dark matter model for the galactic halo. 
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The necessity of introducing dark matter (DM) as the 
main component of galactic matter has become a solid 
fact due to its observational support [l[. Nevertheless, 
the nature of DM is one of the most intriguing mysteries 
in physics. A large variety of particles have been consid- 
ered as the main component of DM in the universe and 
only a few of them are still considered as good prospects 
since they must fulfill several requirements 0. Among 
the survivors, the neutralino and the axion are leading 
candidates [3j. The question we address in the present 
work is the following: If DM is mainly composed of ax- 
ions, what type of astrophysical objects will the axions 
form? 

In order to answer this question, we have solved the 
Einstein-Klein-Gordon (EKG) equations in the semiclas- 
sical limit. The source for the Einstein equations is 
the mean value of the energy-momentum tensor oper- 
ator {T^ w ) of a real, quantized scalar field constructed 
with potential energy density given by 0] 
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It is found that the resulting self-gravitating system, 
the axion star, will have an asteroid-size mass (M ~ 
10~ 16 M Q ) and radius of a few meters. This work im- 
proves our previous attempts [5, 6| by solving numerically 
the EKG system for current allowed values for the axion 
mass m and decay constant f a without the necessity of 
any interpolation. Our findings differ from previous es- 
timates where the effect of the potential energy density 
was either neglected Q or was taken into account with 
a wrongsign in the self-interacting term of the poten- 
tial [1, In the first case, it is known that there is 
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a maximum mass for such self-gravitating system given 
by M max = 0.633 rrip/m, where m is the mass associ- 
ated with the scalar field and m p is the Planck mass. For 
the allowed values of the axion mass, 10~ 5 eV < m < 
10 -3 eV 0, the maximum mass for a self-gravitating 
system with the potential energy density (fTJ) neglected 

lies in the range 1CT 8 MQ < M^° n star < 1CT 5 MQ. 
On the other hand, when the axion is considered to have 
a repulsive self-interacting term, instead of the attrac- 
tive one given by (TT]), the maximum mass will be as big 
as M ~ 10 4 MQ. 



Here we solve the EKG system including a Taylor ex- 
pansion of the potential energy density (TT]) and we ob- 
serve that its inclusion tends to decrease the mass and 
consequently the compactness of the self-gravitating sys- 
tem made of axions. Due to the smallness of the axion 
star's masses they could play the role of scalar field mini- 
MACHOs [llf and they will be the final state of axion 
miniclusters U2] originated in the early universe at the 
QCD epoch (l3j |. Assuming that the axion is the main 
component of DM, the galactic halo will be a collisionless 
ensemble of axion stars and will be indistinguishable to 
the standard CDM scenario since A^-body simulations of 
CDM with ultra-high resolution are insensitive to particle 
mass gr anularity smaller than 10 5 M Q - 10 3 M Q [lllll- 



The paper is organized as follows: in section HI the 
EKG equations for a real, quantized scalar field with a 
Taylor expansion of the potential energy density (TTJ) are 
obtained and are solved for arbitrary values of the axion 
mass m and the decay constant f a - In Sec. |H]we include 
the current values of m and f a and we obtain the mass 
and radius of the axion-stars. We finish section [TT] by 
commenting on some consequences derived in the case 
that the axion-star has the properties calculated here. 
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FIG. 1. Gravitational mass as a function of the central value 
of the scalar field <j(0) for different values of A. Dotted lines 
include only an expansion of the axion potential up to the 
term <1> 4 



I. EINSTEIN-KLEIN-GORDON WITH AN 
AXION-LIKE POTENTIAL 

Since axions are real scalar particles, it is very use- 
ful to remember how self-gravitating systems made of 
spin zero particles are constructed. We will follow the 
method developed by Ruffini and Bonazzola 16]. The 
self-gravitating system arise as a solution of the EKG 
equations: 



Gnu = 87rG < T^v > 



(2) 
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and /i^/ m ' ) are the usual creation (annihilation) opera- 
tors for a particle with angular momentum I, azimuthal 
momentum m and energy E n . These operators sat- 
isfy the usual commutation relations \fi„i m , Mn'i'm'l = 

iVnlm'Vn'l'm'} = „ and If^tlm^n'l'm'} = <W<^'<W- 

With the operator <t, it is now possible to construct the 
energy-momentum tensor operator just by inserting 
the operator <f> into the classical expression for the energy- 
momentum tensor 

Tft = g^d^d^ - ^g^dx^ - 5ZV(ct>) . (6) 

The average (Q\T^ U \Q) is done by considering a state 
\Q) for which all N particles are in the ground state 
(Z = m = 0, n = 1). The ground state satisfies 
/ij" 00 |Q) = 0. This procedure, as pointed out in 
cancels all time dependence on the vacuum expectation 
value (Q\Tfj, u \Q) and, for the case of a free scalar field 
(V(<fr) = ^-$ 2 ), the real quantized scalar field yields the 
same field equations as those obtained by using a classi- 
cal complex scalar field. At this level, the self-gravitating 
system for a real quantized scalar field does not differ 
from a complex classical scalar field, hence, a real quan- 
tized scalar field does not produce the so called "oscilla- 
tons" which are time-dependent. In our case we are 
interested in the axion potential (1). In order to compute 
(T 7 "'), we perform a Taylor expansion of ((TJ, i.e. 
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where □ = {1/ \/—g)d tl [y^ r gg ,J ' u du] and V((j)) is the scalar 
field potential. Here < • • • > denotes an average over the 
ground state of a system of many particles. Its presence 
refers to the fact that we are working in the semi-classical 
limit of the Einstein's equations. We will work with units 
where c = H = 1. In the case of a spherically symmetric, 
static space-time described by 



ds 2 = B{r)dt 2 - A(r)dr 2 - r 2 (sin 2 



dO) , (4) 



it has been shown that such self-gravitating systems are 
fully characterized by the scalar field properties, i.e. the 
mass m of the scalar field and its energy density potential 
V($) [I3,[n!]- The total mass of the resulting object and 
the typical radius depend mainly on these two properties 
of the scalar field. The axion is no exception. To deal 
with the quantum nature of the axion field, we have to 
compute the average {T^ u ) in eq. ©• What is usually 
done is to quantize the scalar field $ — > $ = $ + + $~ 



m ^ 2 1 m 4 1 rrr , 6 



(7) 



We will show that the final results do not depend strongly 
on the number of terms considered in the Taylor expan- 
sion of ([1]). The relevant term that should be consid- 
ered is the self-interacting term $ 4 and the sign it carries 
with itself, which differs from the one considered in Bo- 
son Stars (BS) [8]. With the potential (0, it is possible 
to compute (T{f) by performing the quantization and av- 
eraging procedure previously discussed. The computed 
average of the stress energy tensor is 



(T °) 



(Ti) 



E 2 R 2 


R' 2 


m 2 R 2 


m 2 i? 4 


m 2 R 6 


2B 


~ ~2A 


2 


+ 12/ 2 


144/ a 4 


E 2 R 2 


R' 2 


m 2 R 2 


m 2 R 4 


m 2 R 6 


2B + 


~2A 


2 + 


12/ 2 


144/ a 4 + 


E 2 R 2 


R' 2 


m 2 R 2 


m 2 R 4 


m 2 R 6 


2B 


~2A ~ 


2 + 


12/ 2 


144/ a 4 + 



+ 



(8) 



We have dropped all sub-indexes since, as we have al- 
ready pointed out, we will assume that the axion is in 
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its ground state. We can observe that there is no time 
dependence in Furthermore, there are new numeri- 
cal factors that appear due to the averaging performed 



in [13 . For instance 



($ 4 ) = 2R 4 and ($ b ) = 5R b 



such a way that we cannot recover the original cos($// a ) 
from which we departed. Following a similar procedure 
but now applied to the scalar wave equation ([3]), with a 
potential (J7J and the spherically symmetric metric 
the Einstein-Klein-Gordon system is obtained: 
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, we rewrite the system 
rrn, 



Following standard definitions 

(|9]) in dimensionless variables: x = rm, R = cr/y AitG 
and B = m 2 B/E 2 , and we have found it convenient to 
define the dimensionless self-interaction term 
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By imposing regularity at the origin and flatness at infin- 
ity, system ([9]) is solved using a shooting method. Even 
though the set of equations ([9]) is very similar to the case 
for typical BS, 0, the behavior we found for the fam- 
ily of solutions with zero- nodes is completely different. 
A full set of equilibrium configurations is shown in Fig. 
[TJ where the gravitational mass is plotted for different 
values of ct(0) and A. The equilibrium configurations 
have a maximum mass M max at some er(0) = a c for each 
value of A. But the switch in the potential sign of the 
<fr 4 term produces a significant change in the behavior of 
M max in comparison with standard BS The relation 
Mmax ~ A 1 / 2 is no longer satisfied. Instead of increasing 
M m ax as we increase the value of A, we found a decreas- 
ing M max . This effect is expected since instead of adding 
a repulsive interaction between the particles of the sys- 
tem, the change in the sign due to the cosine-like poten- 
tial ([TJ implies an attractive potential. Thus the total 
number of particles needed to form an equilibrium con- 
figuration that balance the gravitational collapse against 
the quantum pressure is lower than the case of a repul- 
sive potential. One can think that this effect is apparent 
and as soon as the complete potential (JT|) is implemented 
in the EKG system, a different behavior would be seen. 
But the decrease in the mass of the equilibrium configu- 
rations is a robust behavior. The masses of equilibrium 
configurations including up to the fourth power of $ in 
the Taylor series are plotted also in Fig. [TJ illustrating 
this robustness. 

The bigger the value of A, the lower the differences on 
the masses. This is because when A is increased, cr c de- 
creases (and equivalently $(r) where we are interested). 
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FIG. 2. JJgg and compactness for different configurations. 
The dependence of i?gg on A is negligible. For a given value 
of <t(0), the compactness decreases as A increases. 



Then, the true expansion parameter of (JJJ is A$, an it is 
always a small parameter . Another interesting issue is 
that the dependence of the radius R gg (defined as the ra- 
dius where 99% of the total gravitational mass is reached) 
on the value of A is weak, as it is shown in upper panel 
of Fig. [2] Combining the invariance of the radius as A 
increases, with the decrease in the mass, means that the 
self-gravitating system made of a scalar field that has an 
axion-like potential has a lower compactness (2M/Rgg) 
as the self-interaction term increases. This "newtoniza- 
tion" of the system is shown on the lower panel of Fig. 



1. 



II. AXION-STAR 



The previous results where obtained assuming arbi- 
trary values of the mass m of the scalar field associated 
with the axion as well as free values for the decay con- 
stant f a . But the mass of the axion is constrained by as- 
trophysical and cosmological considerations to lie in the 
range 10 -5 eV < m < 10~ 3 eV and the decay constant is 

related to the axion mass by m = 6/ieV (^^f^) HHH 

With these two restriction we have 10 13 < A < 10 17 and 
then, the previous selection of dimensionless variables 
{x,er, A.B}) is now inadequate in order to numerically 
solve the system ([9]). After some frustrated attempts, we 
found that a more suitable set of variables to solve the 
system ([9]) is the following: 
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Since A >> 1, it is natural to think that the result- 
ing axion star will have a small compactness and low 
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FIG. 3. Scalar field and potential for a typical axion star 
with different values of cr(0) and an axion mass of m a = 1.0 x 
l(T 5 eV. 



TABLE I. Masses and .Rgg for the configurations shown in 
Fig.! 



a(0) 


Mass (Kg) 


-R99 (meters) 


density p (Kg/m 3 ) 


5 x 1CT 4 


3.90 x 10 13 


1.83 


6.3 x lO 12 


3 x 1CT 4 


6.48 x 10 13 


2.86 


2.7 x 10 12 


1 x 1CT 4 


1.94 x 10 14 


8.54 


3.1 x 10 11 



mass. So, besides the change in variables (fTT|) . it is con- 
venient to solve for a(x) = 1 — A(x). Solving the system 
© for the new set of variables {x,a(x),a(x), B(x)} we 
obtained typical nodeless configurations for these axion 
stars. Some of them are shown in Fig. [3l where we have 
taken an axion mass m = 10 -5 eV. The total gravita- 
tional mass and the radius -R99, both in physical units, 
for those configurations are shown in Table HI 

A possible scenario emerges with the hypothesis that 
DM is mainly composed by axions. As was already 
pointed out by Kolb and Tkachev [12J, nonlinear effects 
in the evolution of the axion field in the early universe 
may lead to the formation of "axion miniclusters" . These 
miniclusters may relax, due to the collisional 2a — > 2a 
process or by gravitational cooling [21(, and they will 
evolve to a BS. In the present work we have constructed 
those BSs for axion particles by solving the EKG system 
for a real quantized scalar field which is regulated by 
an axion-potential JT]). These self-gravitating systems, 
the axion stars, have very small masses and radii of me- 
ters (Table HI and consequently very low compactnesses. 



The resulting densities are not enough to produce stim- 
ulated decays of the axion to photons since they occur 
when r T mpF e / T /(_Rm*/ n ) > 1 which implies densities 
p > 10 15 Kg/ m 3 for m = 10" 5 eV 0, H3- Typical 
densities for axion stars are shown in table HI 

The galactic halo will be an ensemble of axion stars and 
this picture is not in contradiction with observations since 
the size of axions stars fit into the limits coming from 
microlensing or gravothermal instability 11]. Previous 
studies that construct models for galactic dark matter 
halos out of scalar fields assume that each galactic halo 
is a spherical Bose- Einstein condensate made of an ultra- 
light scalar field (23- 28] . In the present work, the axion 
stars will play the role of the scalar field mini-MACHOs, 
that is, a scenario where the scalar field (the axion) form 
a large number of stable asteroid-sized scalar conden- 
sations which end up clustering into structures similar 
to CDM halos with all their advantages. The stability 
of structures made of scalar fields have been extensively 
studied [29M3TI . and those analysis could be extended to 
our case with potential energy density given by (flj. 

Furthermore, axion stars present similar characteris- 
tics to the recently proposed neutralino stars [32], with 
the advantage that they could survive longer periods of 
time 



If DM is distributed as axion stars, their detection 
will be very difficult. The proposed femptolensing 
to detect axion compact objects [13] is close to its 
lower detectable limit. Another related issue is the low 
number of axion stars around the earth. Assuming for 
instance a Navarro-Frenk- White profile for the galactic 
halo, and a local halo density of 0.3 GeV/cm 3 around 
the Sun, there will be ~ 1 axion star in the volume 
cover between Jupiter and the Sun. Nevertheless, 
another axion properties can shed light on the axion, 
such as the conversion of axions into photons in the 
presence of strong magnetic fields [351 ]. Collisions of 
axion stars with neutron stars [36| will produce flashes 
of light that could be detected by Gamma ray Obser- 
vatories 37). A more detailed analysis of these ideas 
together with a more detailed study of the stability 
of the axion stars could help us to determine if DM 
is mainly composed by scalar field particles as the axions. 
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